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ABSTRACT 

A  computer  simulation  of  the  Brayton  cycle  is  used  to  develop  a 
design  procedure  with  respect  to  minimizing  volume  and  weight  for  the 
counterflow  plate-fin  recuperator.   Based  on  the  Orenda  OT-4  600  horse- 
power gas  turbine,  recuperator  performance  and  dimensional  character- 
istics are  presented  for  an  idealized  equilateral  plate-fin  recuperator 
core  matrix.   The  effects  of  plate  spacings  varying  between  0.1  inches 
and  0.5  inches  on  recuperator  performance  characteristics  are  presented 
over  a  wide  range  of  core  frontal  areas.   Specific  trends  toward  minimum 
volume  and  minimum  weight  plate-fin  recuperator  cores  are  discussed  in 
detail. 
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mass  flow  rate  per  unit  area 

acceleration  of  gravity 

convection  heat  transfer 
coefficient 

horsepower 

specific  enthalpy 

thermal  conductivity 

flow  length 

mass  flow  rate 

pressure 

wetted  perimeter 

heat  transfer  rate 

lower  heating  value  of  fuel 


sq  ft 

sq  ft 

sq  ft 

ft 

ft 

ft 

BTU/(lbm  deg  F) 

BTU/(hr  deg  F) 

lbf 

lbm/(hr  sq  ft) 

(lbm  ft)/(lbf  sq  sec) 

BTU/(hr  sq  ft  deg  F) 

hp 

BTU/lbm 

BTU/(hr  sq  ft  deg  F)/ft 

ft 

lbm/ft 

lbf/sq  ft 

ft 

BTU/hr 

BTU/lbm 


English  Letter  Symbols  (continued) 
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I.   INTRODUCTION 

The  heat  exchanger,  a  device  which  effects  the  transfer  of  heat  from 
one  fluid  to  another,  is  found  whenever  it  is  necessary  to  transfer 
energy  without  mixing  the  fluids  involved.   It  is  felt  by  many  that  our 
future  as  a  power  rich  nation  may  depend  partially  on  the  heat  exchanger 
engineer,  as  the  most  promising  sources  of  power  often  present  the  most 
difficult  applications,   One  such  power  source,  the  gas  turbine,  has  a 
number  of  inherent  advantages  over  the  common  gas  and  diesel  engines; 
among  them  are  compactness,  low  weight,  ease  of  maintenance,  and  multi- 
fuel  capabilities. 

With  the  addition  of  a  heat  exchanger  of  reasonable  size,  the  ther- 
mal efficiency  of  a  gas  turbine  unit  may  be  increased  to  about  the  same 
as,  or  slightly  higher,  than  the  gas  engine  [7].*  The  added  weight  of 
the  exchanger,  however,  could  well  negate  one  of  the  primary  advantages 
of  the  gas  turbine.   Furthermore,  the  relative  rise  in  production  cost 
due  to  the  addition  of  a  heat  exchanger  ought  to  be  considered  in  any 
comparison,  so  as  to  keep  the  gas  turbine  on  a  specific  cost  level  with 
the  diesel  engine. 

Heat  exchangers,  however,  are  inherently  bulky,  heavy,  and  expen- 
sive.  To  overcome  this  problem,  the  concept  of  the  compact  heat  ex- 
changer has  come  into  being.   Not  necessarily  small,  such  an  exchanger 
incorporates  a  heat  transfer  surface  of  high  area  density  (leading  to  a 
high  effectiveness  for  given  weight  and  volume  limitations)  and,  generally, 


* 
Numbers  in  brackets  designate  references  listed  in  Bibliography. 
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may  be  manufactured  at  a  much  lower  cost  per  square  foot  of  heat  trans- 
fer area  than  the  not  so  compact  heat  exchanger  [9]. 

Although  standard  lines  of  compact  heat  exchangers  have  been  adopted 
by  many  manufacturers,  such  exchangers  may  not  be  purchased  as  off-the- 
shelf  items,  as  a  proper  design  is  needed  for  nearly  every  requirement. 
Such  a  design  procedure  usually  considers  the  factors  of  heat  transfer 
requirements,  pressure  drop  characteristics,  size,  and  cost,  with  an 
emphasis  on  the  balance  between  the  gain  due  to  heat  recovery  and  the 
loss  due  to  pressure  drop  contributed  by  the  addition  of  the  exchanger. 

There  exist,  however,  numerous  applications  for  which  this  much 
used  design  philosophy  is  not  sufficient.   Often  the  designer  is  faced 
with  stringent  volume  or  weight  limitations.   The  technique  of  deter- 
mining an  optimized  counterflow  heat  exchanger  with  respect  to  volumetric 
and  weight  considerations  is  developed  in  this  study. 

Fundamental  relationships,  assumptions,  and  design  points  upon 
which  the  investigation  is  based  are  presented  in  Sections  II  through 
IV;  the  computer  program  used  to  simulate  the  regenerated  gas  turbine 
cycle  is  discussed  in  Appendix  B. 
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II.   THE  IDEAL  REGENERATED  BRAYTON  CYCLE 

The  gas  turbine  plant,  to  a  large  extent,  must  have  all  its  units 
designed  with  reference  to  the  cycle  on  which  it  is  to  operate.   The 
designer,  therefore,  must  have  a  knowledge  of  the  particular  cycle  to 
be  considered.   The  most  prevalent  type  of  gas  turbine  in  use  today, 
the  simple  open  cycle,  has  achieved  great  success  due  to  its  versatility 
[6].   With  the  addition  of  regeneration  the  efficiency  of  this  cycle  may 
be  greatly  increased. 

The  Brayton  cycle,  the  basic  thermodynamic  cycle  of  the  gas  turbine, 
(refer  to  Figure  1)  is  characterized  by  taking  air,  as  the  working  fluid, 
from  the  atmosphere  into  a  compressor  where  part  of  the  turbine  power  is 
used  to  raise  the  pressure  before  the  air  enters  the  cold  side  of  a  re- 
cuperator.  The  recuperator  reduces  the  heat  input  to  the  cycle  by  trans- 
ferring heat  from  the  hot  turbine  exhaust  gases  to  the  cooler  air.   A 
further  temperature  increase  occurs  as  fuel  is  burned  directly  in  the 
air  flow  (forming  gaseous  products)  through  the  combustor.   The  gas  then 
enters  a  turbine,  where  it  expands  to  a  lower  pressure  while  releasing 
energy  available  for  work.   Upon  leaving  the  turbine  stage,  the  gas 
enters  the  recuperator  hot  side  and  transfers  much  of  its  otherwise  waste 
heat  to  the  incoming  air,  before  being  exhausted  into  the  atmosphere. 

The  particular  cycle  used  in  this  heat  exchanger  investigation 
is  based  on  the  Orenda  OT-4  gas  turbine,  a  600  horsepower  unit  engine- 
ered for  the  U.  S.  Navy  [10],   The  design  specifications  for  the  Orenda 
unit  are: 
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HP  net  power  output  600  hp 

19.  overall  thermal  efficiency  34% 

PI  compressor  inlet  pressure  2120  lbf/sq  ft 

Tl  compressor  inlet  temperature  560  deg  R 

T4  turbine  inlet  temperature  2200  deg  R 

P6  exchanger  exhaust  pressure  2120  lbf/sq  ft 

P2/P1  compressor  pressure  ratio  4.0 

P3  -  P4        combustor  relative  pressure  drop   0.02 

P3 
O  compressor  efficiency  86% 

n  combustor  efficiency  987o 

if)  turbine  efficiency  86% 

t 

T_  fuel  inlet  temperature  530  deg  R 

q  fuel  lower  heating  value  19000  BTU/lbm 

Other  parameters  which  are  usually  specified  for  the  cycle,  such  as  the 
mass  flow  rate,  regenerator  effectiveness,  and  exchanger  pressure  drops, 
are  not  presented;  they  are,  in  fact,  considered  as  variables  in  the 
optimization  procedure  to  be  discussed. 

A  more  detailed  discussion  of  the  Brayton  cycle  is  outlined  in 
Appendix  A. 
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FIGURE  1    THE  IDEAL  REGENERATED  BRAYTON  CYCLE 
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III.   RECUPERATOR  SURFACE  AND  CORE  GEOMETRY 

For  all  applications  in  which  the  turbine  inlet  temperature  is 
limited  and  a  low  to  moderate  pressure  ratio  is  to  be  used,  the  most 
efficient  method  of  improving  the  economy  of  the  gas  turbine  involves 
the  use  of  a  heat  exchanger,  in  which  hot  turbine  exhaust  gases  give  up 
heat  to  the  relatively  cool  air  from  the  compressor.   Two  basic  heat  ex- 
changer models,  the  recuperator  and  regenerator,  are  in  common  use.   The 
term  recuperator  (which  is  considered  in  this  study)  is  usually  applied 
to  a  continuous  flow  method  of  heat  exchange,  while  the  term  regenerator 
specifies  a  periodic  flow  device,  usually  involving  a  moving  matrix. 

A.    COMPACT  HEAT  EXCHANGERS 

Considerations  involving  heat  transfer  characteristics  of  the  gas 
turbine  recuperative  unit  have  led  to  the  development  of  the  compact 
heat  exchanger.   The  term  compact  does  not  refer  directly  to  a  small 
bulk  or  weight,  but  to  the  heat  transfer  surface  density  of  the  ex- 
changer.  The  most  useful  measure  of  this  density  is  the  compactness, 
fc>  ,  where 

Pj   =  total  transfer  area  of  one  side  of  exchanger 
volume  between  the  plates  of  the  same  side 

Somewhat  arbitrarily,  the  working  limits  of  compactness  for  a  simple 

counterflow  recuperator  may  be  specified  as  1000  sq  ft/cu  ft  as  the 

maximum  working  limit,  and  200  sq  ft/cu  ft  as  the  minimum  limit  before 

leaving  the  realm  of  compactness  [9]. 

An  understanding  of  the  geometric  parameters,  presented  in  Figure 

2,  and  subsequent  relations  most  often  used  is  essential  to  any  compact 

heat  exchanger  design  procedure. 
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COMPACT  COUNTERFLOW  RECUPERATOR  CORE 
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The  following  equations  give  the  relations  between  surface  and 
core  factors  for  one  side  of  the  counterflow  plate-fin  recuperator. 
Subscript  1  refers  to  any  one  side,  and  2  refers  to  the  other  side. 
Factors  without  a  subscript  are  common  to  both  sides. 


a      =   Aj,    -  _At 

P±  V,         b,Lc 


3.1 


VU       - 


LAc 
A 


-ii 


3.2 


d,     = 


_    Ai    _         bx  £>.. 


A/ 


V>x*  bi  +2. a. 


3.3 


,_-     -   Aci  -      Ai  *Vi     — 


■?r 


Sr 


3.4 


Ac,  -  °i  A^ 


Air 


Wl 


L 


3.5 


Ai  =  ^iLA(f 


3.6 


Vs    =  lT(^  )-( Ac,   *   A«)' 


3.7 


Among  the  factors  common  to  both  sides  is  A   ,  the  total  recuperator 
frontal  area.   The  preceding  relations  are  sometimes  written  on  the 
basis  that,  for  the  counterflow  recuperator,  the  hot  and  cold  side 
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frontal  areas  are  defined  separately  such  that  their  sum  equals  A 


fr' 


This  investigation  requires  no  such  separation  of  terms  and  is  based 
on  the  interpretation  of  frontal  area  as  presented  in  Figure  2. 

B.    PLATE  -  FIN  SURFACES 
1.   General  Theory 

Although  other  matrix  forms  are  used,  the  plate-fin  recuperator 
matrix  is  the  most  common  [6],   Not  only  does  this  type  of  structure  lend 
itself  to  light  weight,  extremely  compact  designs,  but  the  plate  spacing 
on  each  side  may  be  easily  optimized  independently,  allowing  for  design 
flexibility. 

An  understanding  of  the  specific  geometric  parameters  of  the 
plate-fin  matrix  (Figure  3)  and  their  relations  is  necessary  to  the 
design  procedure. 


FIGURE  3 
PLATE  -  FIN  GEOMETRY 
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a  plate  thickness  ft 

b  plate  spacing  ft 

o  fin  thickness  ft 

S  fin  spacing  ft 

A  heat  transfer  area  sq  ft 

A  free  flow  area  sq  ft 
c 

P  wetted  perimeter  ft 
w 

L  recuperator  length  ft 

w  fin  thickness  parameter  ft 


Ian© 


A    -  fc,  L 


3.8 


Ac-  -^  SC^*W)  3.9 


Pw=  Zi  ♦  S    -  a(a^-e)-  S    =    s(c-^-&  ♦  i)  3.10 


3.11 


From  these  relationships,  the  values  for  two  very  important  parameters 
may  be  found. 


r\  - 


_   A 


c 


=  (b-w)[2CS7e  +  0] 


h~   ~p^    -  ^o-w;|cvCo^e  '  ±J\  3.12 
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-  &R>  -  ^(cSTe  *  x) 


^=A[(ibi)LJ    '£%  ='8-^*1J 


Since  the  recuperator  matrix  studied  involves  a  fin  angle,  ©  , 
of  sixty  degrees  (resulting  in  an  equilateral  fin  matrix)  a  further 
simplification  may  be  made. 

W  =  2-cT  3.14 

-  (V>-2Q  3  15 

rK g  3.15 

Table  I  represents  values  of  these  parameters  for  the  plate  spacings 
investigated.   Although  the  spacings  greater  than  0.4  inches  fall  out 
of  the  previously  defined  compactness  region  (200  i  b  i  1000) f  they 
are  included  in  the  study  since:  (1)  this  region  is  arbitrary  and  (2) 
recuperator  trends  in  weight  and  volume  are  more  easily  observed  when 
large  and  small  spacings  are  compared. 
2.   Specific  Parameters 

In  order  to  keep  the  investigation  on  a  common  basis,  all  cores 
analyzed  are  assumed  to  be  constructed  of  stainless  steel,  with  con- 
stant properties,  and  uniform  plate  and  fin  thicknesses.   The  constant 
parameters  used  in  the  analysis  are: 

a  plate  thickness  0.0008333  ft 

(0.01  in) 

£  fin  thickness  0.0004166  ft 

(0.005  in) 
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k  thermal  conductivity  of       12  BTU/(hr  sq 

stainless  steel  plates        ft  deg  F)/  ft 
and  fins 

Q  density  of  stainless  484  lbm/cu  ft 

steel  plates  and  fins 
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TABLE  I 

EQUILATERAL  PLATE  -  FIN  SPACINGS 

Plate  Compactness,  P>  Hydraulic 

Spacing, b  Radius, rh 

(In)  (sq  ft/cu  ft)           (ft) 

0.1  720.00  0.0012500 

0.125  576.00  0.0015972 

0.15  1*80.00  0.0019itUit 

0.175  411. U3  0.0022917 

0.2  360.00  0.0026389 

0.25  288.00  0.0033333 

0.3  21*0.00  0.001*0278 

0.35  205.71  0.001*7222 

O.t*  180.00  0.005U167 

0.5  11*1*. 00  0.0068056 


The  hydraulic  radius,  r.,  Is  based 

n 

on  a  fin  thickness,  6*  ,  of  0.005  In. 
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IV.   HEAT  TRANSFER  AND  FLOW -FRICTION  DESIGN  THEORY 

The  design  of  a  recuperator  involves  a  consideration  of  both  the 
heat  transfer  rates  between  the  fluids  and  the  pumping  power  expended  to 
overcome  fluid  friction  in  the  exchanger  [8],   Obviously,  one  of  the  most 
desirable  characteristics  of  a  compact  heat  exchanger  is  that  it  yield  a 
maximum  heat  exchange  for  a  minimum  pressure  drop. 

A.   HEAT  TRANSFER  CHARACTERISTICS 

Elementary  heat  transfer  principles  show  that  the  cool  air  and  hot 
gas  flows  should  be  in  opposite  directions ;  otherwise,  the  two  streams 
would  approach  their  mean  temperature  only.   Subsequently,  many  compact 
heat  exchangers  involve  a  counterflow  design,  with  plate-fin  construc- 
tion for  added  efficiency  [5]. 

An  overall  heat  transfer  rate  equation  which  includes  convective 
and  conductive  heat  mechanisms  responsible  for  the  heat  transfer  from 
the  hot  gases  to  the  cool  air  may  be  written  for  a  differential  element 
as 

dc^  -  U  dA  (Tw-Tc)  4.1 

where  U  is  an  overall  heat  transfer  coefficient,  dA  represents  an 
element  of  heat  transfer  area,  and  dq  is  a  differential  heat  flux. 
A  physical  understanding  of  the  heat  transfer  performance  of  the  heat 
exchanger  may  be  attained  with  the  use  of  equation  4.1  and  a  number  of 
parameters  to  be  defined. 

1.    Fin  Effectiveness 

The  straight  recuperator  fin  of  constant  thickness   o   and 
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length  Ji   is  represented  in  Figure  4.   Such  a  fin  extends  from  the 

base  plate,  at  temperature  T_,  through  the  moving  fluid,  at  temperature 

T       ,  to  a  second  plate.   As  there  is  no  heat  flux  through  the  fin 
cc 

midway  between  the  plates  (due  to  the  symmetry  involved),  a  heat  trans- 
fer analysis  need  only  consider  a  fin  of  effective  length  Si   /2,  insul- 
ated at  the  tip  (Figure  4).   An  energy  balance  over  an  element  of  this 
fin  [2]  results  in  the  equation: 


q  =  k  m  A  0oTanh(m  *■) 


4.2 


k  =   thermal  conductivity  of  fin  material 

=  a  dimensionless  parameter  for  thin  fins 


m 


*   / 


2h 
k£ 


A  £  &(1)   =  cross  section  area  of  fin  per  unit  width 


f 
So 


=  T   -  T 
i0     0° 


h  =  convective  heat  transfer  coefficient  between  the  fin  and 
the  fluid 


FIGURE  4 
STRAIGHT  FIN  OF  UNIFORM  THICKNESS 
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In  the  order  to  express  the  heat  exchanging  capacity  of  an  extended 
surface  it  is  useful  to  define  the  term  fin  effectiveness  as 

n    =  fin  effectiveness 

actual  heat  transferred  from  fin 

^V   "  heat  transferred  if  entire  fin  were  at  temperature  T 

k  m  At  6o   Tanh(m  JL  /2) 
"  t  =        2  L  h  (  A  /2)e0  4.3 

r\l     =   (Tanh(m  Jl/2))  (m  Jl/2)"1  4.4 

2.    Surface  Effectiveness 

Since  the  plate-fin  geometry  involves  an  array  of  fins, 
it  is  useful  to  define  a  total  surface  temperature  effectiveness,  O0  , 
which  gives  a  measure  of  performance  of  the  entire  array.   Defining 
"A"  as  the  total  heat  transfer  area  of  the  combined  fin  and  plate 

surfaces, 

actual   heat    transfer   from  A  4.5 

(•°       heat    transfer   if  A  were   at  T 

XI      =    (A    -  Af)h(TQ    -   T(0>  +     r^      Afh(TQ    -   Tqq    )  ^ 

A     h      (T      -Ttf> 


n    =1-   *±  Cl-ilO 


3.        Overall  Heat  Transfer  Coefficient 

The   overall   conductance  U,    as   defined   previously,   may 
now  be   calculated   as 


Uc  "  L    cak/aO 


fUhw  C Aw/Ac)  k  Ro.chc, 


4.7 


where  U   is  based  on  a  unit  of  cold  side  total  surface  [2].   In  the 
c 
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usual  gas-to-gas  recuperator,  the  second  term  in  equation  4.7  (the  wall 
resistance  component)  may  be  neglected  relative  to  the  fluid  resistances, 


n.o,v\       overall  surface  effectiveness  for  hot  and 
^o.c      cold  sides 

A         average  wall  area 
w 

A,        total  heat  transfer  area  of  hot  and  cold  sides 
h 

A 
c 

h        convection  heat  transfer  coefficients  of 

h         hot  and  cold  sides 
c 

a         plate  thickness 

k         thermal  conductivity  of  plate  material 


4.    Counterflow  Exchanger  Heat  Transfer  Effectiveness 

As  previously  stated,  the  most  common  and  most  effective  flow 
arrangement  for  a  two-fluid  heat  exchanger  is  the  counterflow  design. 
Referring  to  Figure  5,  an  energy  balance  leads  to  a  second  equation  for 
the  amount  of  heat  transferred  between  the  two  fluids. 

dc^  =    CmCp^rx    dTV    =    Uc?k  dTt  4.8 

Integration  over  the  length  of  the  recuperator  leads  to  the  equation 

C^-  Cm.  C^Vl  CTy,  irx  -T^  Qol)  -  U^pV  Ol.oot  -^Aa)   4.9 

It  is  useful  to  define  the  term  "capacity  rate"  as  the  product  of  the 

specific  heat  of  one  fluid  and  the  mass  flow  rate  of  that  fluid,  C, 

referring  to  the  hot  gas,  and  C  referring  to  the  cooler  air. 

c 

c^-    Cn  CTV.wi  -T\,,Oo0    -    C^CTcput    -  Tc/wO  4.10 
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h,in 


c,out 


h,out 


"c.in 


flow  length 

FIGURE  5 
TEMPERATURE  VARIATIONS  IN  COUNTERFLOW  HEAT  EXCHANGER 

The  heat  transfer  performance  of  the  recuperator  may  now  be  deter- 
mined by  the  heat  exchanger  effectiveness,  £  ,  defined  as  the  ratio  of 
the  actual  heat  transfer  to  the  thermodynamic  maximum  heat  transfer 
rate.   The  maximum  heat  transfer  would  be  obtained  if  one  of  the  fluids 

were  to  undergo  a  temperature  change  equal  to  (T,  .   -T    ).  the 

h,in   c,in  ' 

maximum  temperature  difference  in  the  counterflow  exchanger.   Since  the 

heat  lost  by  the  hot  gas  must  be  gained  by  the  cooler  air  (neglecting 

losses),  as  may  be  seen  in  Equations  4.9  and  4.10,  the  fluid  undergoing 

this  maximum  temperature  change  must  have  the  lower  of  the  two  capacity 

rate  values,  C  .  . 
min 


max 


(T 


-  T 


min  h,in    c,in 


) 


4.11 


The  effectiveness  may  now  be  determined  as 

actual  heat  transfer 


£  = 


maximum  heat  transfer 
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£    -       CVT.CTV.tn  -TV.pvA.  J       _         Cc-  CTc  pyA.  -  Tc.|VJV )  4.12 

C«vv.n.  CT^^a  -T£_  La  )     Oaxn.  CTV  ,*un.-T\_ in.) 


From  the  previous  equations,  the  effectiveness  may  be  shown  to  be 
independent  of  the  fluid  temperatures  [2], 

g-     J-  -  -0- 4.13 

i  C_wv.cn.       C  V  -  ?=___£!■  \  C7A 


where  C  .   and  C     are,  respectively,  the  smaller  and  the  larger  of 
min      max     '    r  J  *  & 

the  two  capacity  rate  magnitudes. 

Finally,  defining  the  Number  of  Transfer  Units,  NTU  (a  dimension- 
less  expression  for  the  heat  transfer  size  of  the  recuperator),  as  the 

ratio  of  (U  A)  to  C  .  ,  where  A  refers  to  the  total  surface  area  on 

min 

which  U  is  based,  the  effectiveness  may  be  written  as 

Cr<\A.n. 


t-nv.n.  \ 

•U- £53* )  w™3 


f    -  4  '  -0- 4.14 


A  comparison  of  this  equation  with  effectiveness  relations  for  other 
geometries  will  lead  to  the  conclusion  that  the  counterflow  arrangement 

is  indeed  useful  for  gas  flow  recuperators  (for  which  C  .  /C    *-  1) 

°  mm  max 

when  an  effectiveness  in  the  order  of  80%  is  to  be  obtained. 
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B.    FLOW  -  FRICTION  THEORY 

In  the  design  of  a  gas-to-gas  heat  exchanger,  the  pressure  drop 
characteristics  assume  an  importance  equal  to  that  of  the  heat  trans- 
fer characteristics  [8j.   The  designer  is  generally  interested  in  the 
overall  core  pressure  drop,  including  entrance  and  exit  effects,  rather 
than  the  core  friction  pressure  drop  alone. 


a                                                   b 
I I 

EZZZZZZZZZ2ZZZZZZZZZZZZZ2ZZ2  ZZ2  2  ZZ2  "  A 

Air  Flow  , 

v  ) ,  TZZZZ1  72  I  ZZ2  TE22  TZZZZH  >  >^  TD.  2  t  1  72  "~n 

Hot  Gas  Flow     <- 


177-77-y  f    '  1  j  I   ■'  -  J  f  I  I  J  ZZ2    '  I  U  UOOOOL  I  nj7777A 


FIGURE  6 
RECUPERATOR  CORE  MODEL  FOR  COOL  AIR  FLOW  PRESSURE  DRO] 


L    Core  Friction  and  Flow  Acceleration 


The  core  friction  pressure  loss  term,  due  to  viscous  shear 
and  pressure  force,  is  expressed  in  terms  of  the  dynamic  head   -s  ■—   , 
the  appropriate  friction  factor,  \       ,  and  the  ratio  of  the  length  to 


hydraulic  radius  of  one  side,   — 


AFWWv-  ^  ( Ac  )  C^cf1"  ) 


A   ' 
c 


as 


4.15 
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AP       -_C(A\(_G^\  416 

Friction.    )  {Ac)  ^q  o«^  y 

where  G  is  based  on  the  minimum  free  flow  frontal  area  (  G- p. u»»  -  o^ul?"  tt  ) 

and   P"v  is  the  average  core  fluid  density  (pre-  ^ '1 )• 

Since  density  changes  occur  in  the  core  as  the  fluid  temperature 
changes,  a  term  accounting  for  flow  acceleration  must  be  included  in 
the  pressure  drop  analysis.   The  change  in  momentum  due  to  this  fluid 
acceleration  may  be  expressed  as 

r^  F 

A  P  =  — ■ 4  17 


^R«.u.UTatiw  =  -§-  S~    (    4r^-  -  i^  4.19 


tS-0 


a    pa-  V  pto 

2.    Entrance  and  Exit  Effects 

The  entrance  pressure  drop,  due  to  the  free  flow  area  change 
as  well  as  to  the  irreversible  free  expansion  following  the  abrupt  in- 
let contraction  may  be  written  as  [8]. 

i.  \  u^r 


Px  -  Pa,    =  (  z§)  p,u.vz  U-  <5~l)  +  p,  Kc  2^ 


4.20 


G2 

Px     -    Pcx    -      e^p»      C±-  <^~2     t    Kt)  4.21 
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where  a  constant  density  is  assumed  to  exist  from  point  "1"  to  point 


"a'.   K  ,  the  contraction,  or  entrance,  coefficient  may  be  calculated 
by  methods  explained  in  Ref.  [8]. 

Similarly,  the  exit  pressure  rise  is  due  to  the  change  in  flow 
area  (which  is  identical  to  the  corresponding  entrance  pressure  drop 
term)  and  the  irreversible  free  expansion  and  momentum  changes  follow- 
ing the  abrupt  exit  expansion.   This  pressure  rise  is  written  as 

Z 


u^i* 


Pvs  -  P*_     -  [zq  )  paUL2   Co-2-  l)    t  p%  We  ir^r 


4.22 


P^    -  P2      _     Z^pz 


C  cr-*  -  1  *  We) 


4.23 


where  K  ,  the  exit  coefficient,  is  evaluated  in  Ref.  [8]. 


3.    Overall  Pressure  Drop 

The  combination  of  all  pressure  drop  effects  results  in  the 
equation 


Px  -  P*.  =  2c\px  _(l-o-2  +  W 


4.24 


As   it  has  been  assumed    that     Oa^pi    and       Ot>  ~   p2.      the    relation   for 
the   overall   exchanger  pressure   drop  in   terms   of   specific   volume,   where 


Pl-P 


C-2v, 


2      Z 
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may  be  written  as 


A      Vr 


Vz 


„\  ,    V2.  N  r     *-»■        VTA  vt     .  ?  >v 

(Kc.  +  i-«r-*)  **CvI  -i>    M"Ac"vI    "   vlCl-cr-Wj 

4.25 
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C.    PARAMETRIC  APPROXIMATIONS 

In  the  preceding  heat  transfer  and  flow  friction  analysis,  a  number 
of  important  engineering  parameters  and  fluid  properties  are  used  in 
fundamental  relationships.   Usually  presented  only  in  tabular  or  graphi- 
cal form,  these  parameters  may  not  be  directly  obtained  in  a  computer 
analysis.   To  overcome  this  problem,  functional  relationships  approxi- 
mating the  tabulated  values  must  be  made  a  part  of  the  computer  program. 

For  the  program  [1],  such  relationships  were  formulated  by  making  a 
"least  squares  fit"  of  the  curves  representing  tabulated  data  of  the 
required  parameters.   The  resulting  functional  relationships,  of  the 
form 

PC.O  -  A  i   +  A2  X.   *  Ai  x.2   + v  A\.  X'"1       4.26 

where:   P(x)  =  desired  parameter 

x  =  variable  on  which  the 
parameter  is  based 

A.  =  polynomial  coefficients 
resulting  from  least 
squares  fit 

i-1  =  order  of  resulting 
polynomial 

are  used  extensively  in  the  computer  analysis.   Table  II  specifies  the 
parameters  for  which  such  polynomial  approximations  are  used,  the 
variable  on  which  each  is  based,  and  the  order  of  the  polynomial. 

Table  II  indicates  that  two  parameters  for  which  functional  relation- 
ships must  be  obtained  are  dependent  on  the  Reynolds  number,  Oul  (4  r.)/jjk_ 
The  Fanning  friction  factor,  \     ,  is  used  in  equation  4.15  and  subsequent 
relations  to  obtain  the  pressure  drop  due  to  core  friction  through  the 
recuperator.   The  Colburn-j  factor,  j,  has  been  obtained  experimentally 
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for  a  number  of  recuperator  geometries  [8]  and  is  used  to  obtain  the 
convection  heat  transfer  coefficient,  h,  as  follows: 
This  heat  transfer  parameter  is  defined  as 

2/3 
j  =  St  Pr  4.27 

Since  j  and  Pr  are  determined  from  functional  relationships,  the 
Stanton  number  may  be  calculated  as 

-2/3 
St  =  j  Pr  '  4.28 

Thus,  the  convection  heat  transfer  coefficient,  used  in  equation  4.5 
and  subsequent  relations,  is  obtained  as 

h  =  St  c   p   u  4.29 

P  J 

Experimental  investigations  of  \    and  j  provide  curves  similar 
to  those  presented  in  Figure  7.   As  indicated,  the  laminar  flow  regime 
(Re  <      1000)  produces  linear  relationships  for  the  ideal  equilateral 
triangular  passage  such  that 

f  -  13.33/Re  4.30 

j  =  2.64/Re  4.31 

The  turbulent  region  (Re  >  5000)  data  may  be  closely  approximated  by 
the  linear  relationship  resulting  from  circular  tube  flow  as 

f  =  0.046  Re'°'2  4.32 

j  =  0.023  Re~°'2  4.33 
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FIGURE  7 


EXPERIMENTAL  CURVES  FOR  -f  AND   j 
The  transition  region  (1000  £  Re  £  5000)  is  also  based  on  experi- 
mental data  and  must  be  approximated  by  a  polynomial  equation  for 
computer  analysis.   The  polynomial  coefficients  for  the  resulting 
equations  approximating  this  transition  region,  as  well  as  those  for 
the  other  parameters  in  Table  II,  are  presented  in  Appendix  B. 
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TABLE  II 

FUNCTIONAL  RELATIONSHIPS 

Parameter         Function  Of         Order  Of 

Polynomial 

ij.     Fanning  friction    Reynolds  number  10 
factor 

j     Colburn-j  factor    Reynolds  number  10 

Pr     Prandtl  number     temperature  9 

ju.     viscosity          temperature  9 

Cp    specific  heat      temperature  5 

i      enthalpy           temperature  6 

Prej   relative  pressure   temperature  8 

T     temperature        relative  pressure  8 

T     temperature        enthalpy  9 
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V.   RESULTS 

A.    OBJECTIVES 

For  many  years  engineers  have  designed  heat  exchangers  by  first 
assuming  core  dimensions  and  then  calculating  mass  flow  rates,  heat 
transfer  rates,  and  pressure  drops  which  would  occur  if  fluids  were 
flowing  through  the  assumed  core.   Extensive  investigations  of  various 
heat  transfer  surfaces  offer  the  designer  enough  latitude  in  surface 
geometry  to  successfully  meet  the  heat  transfer  and  pumping  power  re- 
quirements for  a  specific  application. 

Texts  frequently  refer  to  the  possibility  of  optimizing  the  inter- 
nal core  dimensions  of  a  heat  exchanger  once  a  particular  surface  is 
specified,  yet  very  little  research  has  been  published  in  this  area. 
Such  a  technique  forms  the  basis  for  this  study. 

The  computer  program  [1]  used  in  the  study  simulates  the  open  gas 
turbine  regenerated  cycle  under  constant  load  conditions.   All  adjust- 
ments in  internal  exchanger  geometry  are  reflected  directly  in  terms  of 
the  cycle  thermal  efficiency  and  the  overall  core  size  necessary  to 
obtain  that  efficiency.   As  is  outlined  in  Appendix  C,  any  computer 
simulation  of  the  gas  turbine  cycle  requires  set  tolerance  limits  on 
various  cycle  design  points  which  must  be  met.   Early  investigations  of 
general  recuperator  trends  in  the  study  involved  a  tolerance  limit  of 
3.3%  for  the  desired  thermal  efficiency  of  0.34  and  a  limit  of  0.2%  for 
the  desired  horsepower  output  of  600.   More  detailed  analysis  (as  pre- 
sented in  Figures  15  through  22  and  Tables  V  through  VII)  required 
tighter  tolerances  of  1.45%  for  the  thermal  efficiency  and  0.09%  for 
the  horsepower.   The  use  of  two  tolerance  levels  slightly  affects  out- 
put data  for  the  same  core,  as  may  be  noticed. 
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B.    GENERAL  RECUPERATOR  CHARACTERISTICS 

Trends  in  external  dimensions  of  the  plate-fin  recuperator  due  to 
changes  in  internal  spacings  may  be  noted  by  investigating  a  number  of 
different  cores.   A  sufficient  number  of  cores  may  be  obtained  by  pair- 
ing the  plate  spacings  presented  in  Table  I.   For  convenience  a  parti- 
cular core  will  be  referred  to  in  terms  of  the  ordered  pair  (b  ,b,  )  where 

c   n 

b   represents  the  cold  side  spacing  and  b   represents  the  hot  side  spac- 
ing of  the  core. 

Due  to  addition  of  fuel  in  the  combustor,  the  mass  flow  rate  through 
the  recuperator  hot  side  is  larger  than  the  cold  side  flow  rate  by  a 

factor  of  the  fuel  to  air  ratio,  r,;  .   The  hot  side  mass  flow  rate  must 

'   f /a 

be  accommodated,  then,  by  an  increase  in  fluid  density,  flow  area,  or 
fluid  velocity  (m  -  p  A  v) .   Since  fluid  density  decreases  as  tempera- 
ture increases  the  engineer  is  left  with  only  two  alternatives. 

A  velocity  increase  will  not  only  accommodate  the  increased  mass 
flow  rate,  but  it  will  also  increase  the  heat  transfer  rate  per  unit  of 
surface  area  by  a  factor  something  less  than  the  first  power  of  the 
velocity.   The  friction  power  expenditure,  however,  will  increase  by  a 
factor  between  the  square  and  the  cube  of  the  velocity  [8].   It  is  this 
effect  which  serves  as  the  basis  for  the  design  restriction  that  the  hot 
side  plate  spacing  never  be  less  than  the  cold  side  spacing. 

General  recuperator  volume  trends  may  be  noted  in  Figure  8  for  five 
specific  cores.   Each  point  on  every  curve  represents  the  external  core 
dimensions  which  will  satisfy  the  specified  design  points.   As  might  be 
expected,  from  the  definition  of  the  compact  heat  exchanger,  larger  cores 
are  required  as  one  moves  from  the  more  compact  to  the  larger  spacings; 
for  any  one  frontal  area,  more  length,  and  thus  more  heat  transfer  area, 


40 


is  required  as  spacings  increase.   Excluding  the  extremely  compact  core, 
(0.1,  0.1),  points  of  minimum  volume  for  each  core  fall  below  a  frontal 
area  of  five  square  feet.   Furthermore,  every  core  requires  a  frontal 
area  of  at  least  two  square  feet  before  the  specified  cycle  design  points 
are  met. 

A  more  detailed  view  is  presented  in  Table  III,  in  which  the  minimum 
volume  of  the  (0.1,0.1)  core  is  seen  to  fall  within  a  frontal  area  range 
between  30  and  40  square  feet.   Definitely,  the  minimum  volume  of  such  a 
core  does  not  present  the  optimum  package.   Similar  characteristics  were 
noted  only  for  two  other  extremely  compact  cores,  (0.1,  0.125)  and  (0.125, 
0.125). 

An  increase  in  length  over  a  region  of  increasing  frontal  area  is 
noticeable  in  all  cores  except  those  having  a  cold  side  spacing  of  0.1 
inches.   Occurring  at  a  frontal  area  greater  than  that  for  the  minimum 
volume  core,  this  effect  becomes  more  apparent  as  core  spacings  increase. 
Table  III  indicates  that  the  mass  flow  rates  for  these  cores  decrease  as 
frontal  area  increases.   Not  shown  is  the  fact  that  the  recuperator  hot 
and  cold  side  Reynolds  numbers  also  decrease  with  the  mass  flow  rate. 
The  first  knee  of  each  curve  (that  point  at  which  L  begins  to  increase 
with  Af  )  occurs  when  the  two  Reynolds  numbers  approach  a  magnitude  of 
2000.   Reference  to  Figure  7  indicates  a  decrease  in  the  Colburn-j  fac- 
tor with  decreasing  Re  at  this  point.   Evidently,  the  resulting  decrease 
in  the  convection  heat  transfer  coefficient,  h  (which  is  proportional  to 
j),  is  accompanied  by  an  increase  in  heat  transfer  area  (as  L  increases 
with  A.  )  in  order  that  the  desired  cycle  thermal  efficiency  be  main- 
tained.  This  effect  continues  as  j  decreases  until  the  second  knee  at 
which  the  two  Reynolds  numbers  have  approached  a  magnitude  near  1000. 
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At  this  point,  since   j  begins  to  increase  with  increasing  frontal  area, 
less  heat  transfer  area  is  required,  and  L  begins  to  decrease. 

The  effect  is  not  noticeable  for  those  cores  of  low  plate  spacings, 
as  both  the  hot  and  cold  side  Reynolds  numbers  of  such  cores  fall  in  the 
laminar  region  for  all  frontal  areas.   Reynolds  numbers  for  cores  of  low 
cold  side  spacing  and  large  hot  side  spacing  (i.e.,  (0.1,  0.5))  however, 
tend  to  have  a  maximum  Reynolds  number  near  2000  at  low  frontal  areas, 
followed  by  a  rapid  drop  as  frontal  area  increases.   The  slight  length 
increase  required  by  such  cores  is  not  noticeable  due  to  the  magnitude 
of  the  length  increments  used  at  this  point  of  the  investigation. 

Relative  pressure  drops  on  both  sides  of  all  cores  decrease  to  zero 

as  frontal  areas  increase.   In  the  region  of  minimum  volume  (A,_  <■  5  sq 

li- 
ft), the  cold  side  relative  pressure  drop  appears  to  be  negligible  for 

all  cores  except  that  which  combines  a  small  cold  side  spacing  and  a 

large  hot  side  spacing.   The  hot  side  relative  pressure  drop,  on  the 

other  hand,  appears  to  be  high  only  for  those  cores  in  which  the  two 

plate  spacings  are  of  equal  magnitude. 

C.    MINIMUM  VOLUME  AND  WEIGHT  CORES 

Figure  9  indicates  dimensional  requirements  for  various  cores  in  the 
minimum  volume  region.   As  is  to  be  expected,  required  lengths,  and  there- 
fore volumes,  increase  as  core  spacings  increase.   The  effect  of  forming 
cores  of  unequal  spacings  as  compared  to  equally  spaced  cores  is  presented 
in  Figure  10. 

To  this  point  a  trend  has  developed  in  which  smaller  spacings  re- 
quire smaller  cores.   Based  on  this  trend,  one  might  incorrectly  reason 
that  the  minimum  volume  for  a  given  hot  side  spacing  may  be  obtained 
simply  by  decreasing  the  cold  side  spacing  to  a  minimum  allowable  magni- 
tude.  In  fact,  for  the  cores  of  0.25  inch  hot  side  spacing  this  appears 
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to  be  the  case  as  b   is  reduced  from  0.25  inches  to  0.1  inches.   How- 

c 

ever,  the  trend  presented  by  those  cores  for  which  b,  is  0.5  inches  as 

n 

b  decreases  should  be  noted.   A  decrease  from  a  cold  side  spacing  of 
0.5  inches  to  0.25  inches  reduces  core  volume;  a  further  decrease  to  a 
b   of  0.1  inches  produces  a  core  of  extremely  large  volume.   This  core, 
which  combines  the  large  length  and  frontal  area  requirements  of  the  two 
range  limiting  spacings,  results  in  the  largest  minimum  volume  point  re- 
corded of  all  cores  investigated. 

Similar  analysis  shows  the  effect,  presented  in  Figures  11  through 
13,  of  increasing  the  hot  side  spacing  for  a  constant  b  .   For  all  cores 
presented,  the  b,  increase  results  in  increased  frontal  area  and  length 
requirements;  thus,  the  minimum  volume  for  each  series  of  constant  cold 
side  spacing  cores  increases  as  the  hot  side  spacing  is  increased. 

Recuperator  dimensional  and  performance  characteristics  at  minimum 
volume  points  are  presented  in  Table  IV  ;  each  set  of  cores  listed  con- 
sists of  a  constant  b  with  increasing  b.  .   The  optimum  core  with  re- 

c  h 

spect  to  volume  in  each  set  occurs  at  or  near  that  core  configuration 
for  which  b   is  twice  b  .   The  optimum  core  with  respect  to  weight  for 
each  set  consists  of  the  same,  or  a  slightly  higher,  b,  as  the  optimum 
volume  core.   Such  a  trend  might  lead  one  to  believe  a  rule  of  thumb 
often  used  by  manufacturers  that  the  core  hot  side  spacing  be  twice  the 
cold  side  spacing.   Analysis  will  show  this  to  be  an  inadequate  design 
procedure.   Further  investigation  in  this  area,  however,  should  be  pre- 
ceded by  mention  of  the  performance  characteristics  of  the  minimum  core 
points  as  tabulated  in  Table  IV. 

For  any  one  set  of  cores,  the  cold  side  relative  pressure  drop  in- 
creases with  increased  b,  ,  while  the  hot  side  relative  pressure  drop 

h 
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decreases.   This  effect  is  presented  graphically  for  one  set  of  cores  in 

Figure  14.   The  large  magnitude  of  the  hot  side  relative  pressure  drop 

for  cores  of  b,  £   b   is  to  be  expected  since  the  hot  side  mass  flow  rate 
h      c 

must  be  accommodated  by  an  increase  in  either  flow  area  or  velocity  over 

that  of  the  cold  side  flow.   Since  there  is  little  or  no  area  increase 

(b,  X    b  ),  velocities  on  the  hot  side  of  such  cores  must  be  higher, 
he 

yielding  relatively  large  pressure  drops.  As  mass  flow  rates  generally 

fall  and  flow  areas  increase  as  b,  is  increased  for  each  set  of  cores, 

n 

the  hot  side  relative  pressure  drop  falls.   Evidently,  the  cold  side 
relative  pressure  drop  is  affected  more  by  the  general  length  increase, 

as  b,  increases  for  each  set,  than  by  the  decreasing  trend  in  the  mass 

h 

flow  rate.   Thus,  the  cold  side  relative  pressure  drop  rises  with  b   in 
each  set  of  cores. 

D.    OPTIMUM  CORES 

Concentrating  once  again  on  dimensional  aspects,  Figures  15  through 

19  represent  the  effect  of  increasing  b.  for  a  constant  b   in  the  region 

n  c 

of  optimum  core  spacings  with  respect  to  volume  and  weight.   It  is  ap- 
parent from  these  plots  that  an  optimum  design  procedure  is  not  that  for 

which  b,  ■  2b  .   Tables  V  and  VI  present  the  characteristics  of  cores 
n     c 

optimized  with  respect  to  volume  and  weight;  Table  VII,  which  consists 

of  performance  data  for  cores  of  b,  =  2b  ,  is  presented  for  comparison. 

n     c 

A  comparison  of  Tables  V  and  VI  indicates  that  the  minimum  volume  core 

for  any  one  b   involves  a  higher  mass  flow  rate  and  less  length  than  the 
c 

corresponding  minimum  weight  core.  This  results  in  a  lower  cold  side 
and  higher  hot  side  relative  pressure  drop  through  the  minimum  volume 
core  as  compared  to  the  minimum  weight  core  for  any  set  b  . 
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Further  investigation  of  the  Tables  indicates  that  the  minimum  weight 

core  for  any  b   requires  a  larger  b   than  that  for  the  minimum  volume 

core.   This  effect  (noticed  earlier  in  Table  IV)  is  to  be  expected,  as 

the  minimum  weight  core  represents  a  compromise  between  a  low  core  volume 

and  a  low  overall  core  density  (which  decreases  as  spacings  increase). 

Differences  between  the  two  sets  of  optimized  cores  and  those  cores  for 

which  b,  =  2b  are  illustrated  in  Figure  20. 
h     c 

Allowing  a  tolerance  limit  of  0.004  inches  either  side  of  the  optimum 
hot  side  spacing  for  any  given  cold  side  spacing  as  presented  in  Tables 
V  and  VI,  smooth  optimum  spacing  curves  may  be  drawn,  as  shown.   In  the 
range  of  cold  side  spacings  between  0.1  and  0.225  inches  a  linear  re- 
lationship closely  approximates  each  of  the  two  optimum  spacing  curves, 
as 

bu  -  0.86  b   +  0.12  (0.1  ±.   b   1  0.225)   5.1 

he  c 

for  optimum  volume  spacings  and 

bu  =  1.16  b   +  0.12  (0.11  b   £  0.225)    5.2 

he  c 

for  optimum  weight  spacings. 

Outside  this  range  of  cold  side  spacings  the  optimum  core  plots 
tend  to  level  off,  thus  tending  to  keep  the  cores  within  the  limits  of 
compactness  of  the  compact  heat  exchanger  region. 

Since  the  engineer  is  usually  at  liberty  to  choose  both  plate  spac- 
ings, he  must  have  an  idea  of  the  relative  volumes  and  weights  between 
the  optimum  cores  which  he  may  use.   Figure  21  presents  volume  character- 
istics of  optimized  cores  from  Table  V  with  respect  to  the  minimum  volume 
core  (0.125,  0.225).   As  may  be  seen,  if  the  engineer  expects  to  use  a 
cold  side  spacing  equal  to  or  greater  than  0.175  inches,  he  must  accept 
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the  fact  that  his  optimum  core  will  have  a  volume  at  least  20%  greater 
than  the  minimum  available „   Similarly,  Figure  22  represents  the  weight 
characteristics  of  cores  from  Table  VI  with  respect  to  the  minimum  weight 
core  of  (0.175,  0,332) .  Although  the  variance  in  weight  is  not  as  great 
as  that  for  volume,  on  a  percentage  basis,  it  may  be  seen  that  the  engineer 
must  expect  weights  at  least  5%  greater  than  the  minimum  available  if  he 
wishes  to  use  cold  side  spacmgs  other  than  those  between  0=,125  and  0.235 
inches.   The  very  large  increase  in  relative  weight  as  extremely  compact 
cores  (b  i:  0.1  inches)  are  optimized  should  also  be  noted. 

For  the  design  points  of  this  study,  if  the  engineer  were  free  to 
optimize  any  non-dimensional  parameter  (as  mass  flow  rate  or  relative  pres- 
sure drops)  while  keeping  within  10%  of  the  minimum  volume  possible,  refer- 
ence to  Figure  21  and  the  use  of  equation  5-1  would  indicate  a  freedom  of 

b   choice  between  0*1  and  0,15  inches  and  b,  choice  between  0.206  and  0.28 
c  h 

inches.   A  similar  design  procedure  with  a  limit  of  57o  over  the  minimum 

weight  core  (using  Figure  22  and  equation  5-2)  would  yield  a  b   range  be- 

c 

tween  0.13  and  0,225  inches  and  a  b,  range  between  0.27  and  0.38  inches. 

h 

A  reference  back  to  Table  IV  will  indicate  that  the  minimum  weight  and 
minimum  volume  cores  for  all  cores  investigated  do  fall  within  the  indi- 
cated ranges  of  spacings.   If  both  volume  and  weight  of  the  core  are 
limited  within  the  specified  107*  and  57e9  the  range  of  allowable  cold  side 
spacings  is  further  reduced  to  between  0.13  and  0.15  inches,  and  the  allow- 
able hot  side  range  to  0,27  to  0.28  inches. 

Although  the  design  limitations  outlined  are  applicable  only  to  the 
specific  core  geometry  studied,  the  procedure  may  be  applied  to  all  simi- 
lar geometries.   As  indicated,  when  core  and  weight  tolerances  are  speci- 
fied, the  engineer  does  have  some  freedom  in  choosing  internal  core 


46 


dimensions  based  on  other  requirements.  In  all  applications  for  which 
the  allowable  volume  or  weight  of  the  gas  turbine  is  at  a  premium,  it 
seems  that  the  engineer  would  be  wise  to  follow  a  design  procedure  for 
the  heat  exchanger,  as  outlined,  which  optimizes  the  limited  quantity, 
before  moving  on  to  further  core  optimization  with  respect  to  perform- 
ance. 
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APPENDIX  A 

THE  REGENERATED  BRAYTON  CYCLE 

A.    THERMODYNAMICS 

Ideally,  the  basic  open  cycle  consists  of  the  four  processes  of 
isentropic  compression,  constant  pressure  heat  addition,  isentropic 
expansion,  and  constant  pressure  heat  rejection  [4].   The  ideal  regener- 
ated cycle  differs  only  in  that  the  heat  addition  process  is  accomplished 
through  the  use  of  a  recuperator  as  well  as  the  standard  combustor.   With 
these  points  in  mind,  the  ideal  cycle  (refer  to  Figure  23)  may  be  des- 
cribed using  the  equalities; 

A.l 

A. 2 
A. 3 

A. 4 

The  most  obvious  differences  between  the  actual  and  ideal  cycles, 
from  the  thermodynamic  viewpoint,  are  the  increases  of  entropy  and  drops 
in  pressure  due  to  actual  working  conditions.   These  differences  may  be 
best  explained  by  taking  the  cycle  in  its  respective  stages. 

1.    Compression 

State  point  "2n  lies  to  the  right,  at  an  increased  entropy,  of 
ideal  point  9,2a8'  due  to  the  irreversibility  of  the  compression  process; 
both  points,  however,  lie  on  a  line  of  constant  pressure,  since  the 
compressor  outlet  pressure  is  specified  by  working  conditions.   The 
resulting  equality  is 


P2a 

P3a  = 

P4a 

P5a 

= 

P6a 

T2a 

= 

6a 

T3a 

= 

T_ 
5a 
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k.  s 
3  > 

*■>  a 

«D  — 

I.  (TJ 

o>  x: 
aw 
E  c 


2a 


1,1a 


>  5 


Entropy, s 


("a"  Indicates  ideal 
state  points) 


Recuperator 


FIGURE  23    THE  ACTUAL  AND  IDEAL  REGENERATED  BRAYTON  CYCLES 
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P2a  =  P2  A-5 

2.  Heat  Addition 

Two  changes  are  apparent  when  one  compares  cycle  conditions  at 

the  recuperator  exit  point.   First,  the  actual  pressure  level  is  below 

that  of  the  ideal,  due  to  the  pressure  drop  as  the  air  flows  through 

the  recuperator. 

P,  <  P-  A. 6 

3    3a 

Second,  the  temperature  of  the  actual  case  is  below  that  for  the  ideal, 
since  less  heat  is  transferred  than  for  the  reversible  process. 

T,  <  T_  A. 7 

3  3a 

In  order  to  raise  the  temperature  further,  the  air  flows  through  a 
combustor.   The  inlet  temperature  of  the  turbine  is  specified  by  its 
material  properties, 

T.  -  T,  '•       A.  8 

4a    4 

However,  point  "4"  does  lie  at  an  entropy  state  greater  than  point  "4a" 

due  to  the  irreversibilities  involved  between  the  compressor  and  the 

turbine.   Furthermore,  an  additional  pressure  drop  occurs  as  the  air 

flows  through  the  combustor. 

P.  <    P,  A. 9 

4  4a 

3.  Expansion 

An  additional  entropy  increase  results  from  the  expansion 
process  in  the  turbine,  moving  point  "5"  to  the  right  of  the  ideal  state 
point.   As  the  turbine  outlet  pressure  is  not  specified,  and  a  further 
pressure  drop  is  expected  as  the  gases  pass  through  the  recuperator, 
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the  pressure  at  point  "5"  is  necessarily  greater  than  that  for  the 
ideal  cycle. 

P.  v   P.  A.  10 

5  5a 

4.    Heat  Rejection 

As  indicated,  the  pressure  drop  experienced  by  the  hot  gases 
in  passing  through  the  recuperator  results  in  the  pressure  at  point  "6" 
being  equal  to  the  ideal  pressure,  since  this  point  is  considered  to  be 
the  exhaust  to  atmospheric  conditions. 

P,   -  P,  =   14.696  psi  A.  11 

6a     6  r 

Finally,  since  less  heat  is  transferred  from  the  hot  gases  to  the  in- 
coming cold  air  as  is  ideally  possible,  the  temperature  at  point  "6" 
is  greater  than  that  for  the  ideal  case. 

T.  7  T,  A. 12 

6  6a 

B,    CYCLE  SYSTEMS 

An  understanding  of  the  efficiencies  of  each  system  is  necessary 
if  one  is  to  work  with  the  Brayton  cycle.   This  may  best  be  accomplished 
by  looking  at  each  system  process  individually.   It  must  be  noted  that 
the  assumption  of  steady  flow  and  no  potential  energy  change  of  posi- 
tion is  made  for  the  cycle.   More  detailed  information  may  be  obtained 
in  Refs.  [3]  and  [5] . 

1.    Compressor 

r^c  =  compressor  efficiency 

P   _  work  of  isentropic  compression  =   la   2a 
actual  internal  work    i-  -i 

A. 13 
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2.    Recuperator 

C  =  recuperator  effectiveness 


4    = 


actual   heat   transfer 

maximum   possible   heat    transfer 


£     =    C^<V)h  CT\\n.  -^.oolj    -(mc^cCT^d-lA    14 


A  more  detailed  explanation  of  recuperator  effectiveness  is 
presented  in  Section  IV. 

3.    Combustor 

0.    =  combustor  efficiency 

p.     _        actual  energy  received  by  gases 

(uji~i&   energy  released  during  complete  combustion 

<l««»-  ^^f/^U-ia  -g^iy  A#15 

This  definition  may  be  more  easily  understood  when  an  energy  bal- 
ance based  on  one  pound  of  air  and  £  $./<x  pounds  of  fuel  per  pound  of 
air  entering  the  combustor  is  made. 


r;/o.H    +  *W  <U   +  ^  «un*,ou    ^i-fu^*)r^  <U 


A. 16 
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Hi- ^co^b)  r.$./0_  q  $_ 


>  <*  i 4  V  ^) 


FIGURE  24 
COMBUSTOR  ENERGY  DIAGRAM 

of  the  three  entering  energy  terms,  it  can  readily  be  seen  that  i   (the 
enthalpy  of  the  incoming  air)  and  r  .  if   (the  fuel  to  air  ratio  times 
the  enthalpy  of  the  fuel)  represent  molecular  internal  energy  and  flow 
work.   The  third  term  represents  the  chemical  energy  stored  in  the  fuel, 
rf/   q. ;  since  complete  condensation  of  water  vapor  in  the  air  flow  will 
not  occur,  the  lower  heating  value  of  the  fuel  (which  disregards  the 
latent  heat  of  condensation),  q. ,  is  usually  used  in  gas  turbine  anal- 
ysis [3] . 

The  fuel  to  air  ratio,  a  very  important  parameter  in  the  gas  tur- 
bine analysis,  may  now  be  shown  to  be 


A. 17 


4.    Turbine 


It 


turbine  efficiency 

actual  internal  work 

work  of  isentropic  expansion 

1  ^(X.      ~     I^jOl. 


A.  18 
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5.   Additional  Calculations 

Three  final  calculations  are  frequently  useful  in  a  complete  anal' 
ysis. 

Overall  Thermal  Efficiency 

Oa.   -  overall  thermal  efficiency 

O    -    net  work 
"^    heat  supplied 


A.  19 


Horsepower 

Using  the  conversion  factor  2544  BTU  =  one  horsepower- 
hour,  the  equation  for  the  horsepower  output  may  be  seen  to  be 

HP  =  Cfl  tViXT U<x) Iq,  1 ) Cj?l^  A. 20 

2544 

Specific  Fuel  Consumption 

fuel  burned 


SFC  = 


net  work 


SFC  =  tt*/^*°  A. 21 

HP 

2544 
SFC  =  A. 22 

qi 
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APPENDIX  B 
COMPUTER  PROGRAM  ESSENTIAL  FEATURES 

A.    GENERAL  ASPECTS 

The  computer  program  referred  to  in  Ref.  [1]  has  been  greatly  ex- 
panded by  Professor  P.  F.  Pucci  since  its  original  use,  and  may  now  be 
used  to  investigate  parallel  and  crossflow,  as  well  as  counterflow, 
recuperator  characteristics.   This  counterflow  analysis  is  based  on 
design  parameters  for  the  Orenda  OT-4  gas  turbine  only;  the  program, 
however,  is  applicable  to  any  gas  turbine  recuperator  analysis,  as  the 
desired  design  points  are  introduced  as  input  data.   A  basic  under- 
standing of  the  program  logic  may  be  obtained  by  referring  to  the  simpli- 
fied flow  chart,  Figure  25;  a  sample  program  output  is  presented  in 
Figure  28. 

The  program  is  written  on  the  basis  that  changes  in  recuperator  core 
geometry  (plate  spacings  and  fin  configurations)  and  dimensions  (frontal 
area  and  length)  are  reflected  directly  in  the  overall  cycle  thermal 
efficiency,  rather  than  in  the  often  used  design  variables  of  heat  ex- 
changer pressure  drop  and  heat  transfer  characteristics.   Although  the 
program  may  be  used  to  obtain  data  for  a  number  of  conditions  (such  as 
the  computation  of:   (1)  the  maximum  cycle  thermal  efficiency  given  a 
recuperator  of  specified  volume,  (2)  recuperator  dimensions  necessary  to 
obtain  a  specified  thermal  efficiency  and,  (3)  general  outputs  for  a 
simple  open  cycle,  to  mention  only  a  few)  only  two  output  types  were  used 
extensively  in  the  analysis. 
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The  first,  Output  Type  III,  computes  the  required  recuperator 
length  for  a  specified  frontal  area  and  internal  geometry;  furthermore, 
the  horsepower  and  cycle  thermal  efficiency  are  specified,  within  pre- 
determined tolerance  levels.   This  output  is  extremely  useful  in  provid- 
ing an  overall  picture  of  the  desired  recuperator's  behavior  over  a  wide 
range  of  frontal  areas.   Information  on  the  effect  of  changes  in  internal 
core  dimensions,  may  be  obtained  with  the  use  of  Output  Type  VIII.   This 
routine,  similar  to  Output  Type  III,  computes  the  required  length  of  a 
recuperator  over  a  range  of  frontal  areas.   In  this  procedure,  however, 
the  cold  side  plate  spacing  remains  constant,  as  before,  whereas  the  hot 
side  plate  spacing  is  incremented  over  a  desired  range.   Flow  charts  for 
these  routines  are  presented  in  Figures  26  and  27. 

B.    PARAMETRIC  APPROXIMATIONS 

As  explained  in  Section  IV,  the  computer  program  makes  use  of  poly- 
nomial approximations  to  obtain  a  number  of  fundamental  parameters.   The 
polynomial  coefficients  used  in  this  analysis  will  be  presented,  should 
a  need  for  them  arise  in  the  future;  unless  otherwise  indicated,  the 
polynomial  is  of  the  form 

P(x)  =  A  +  Ax  +  Ax  +...............+  A. x  B.l 

where:  P(x)  =  desired  parameter 

x  =  variable  on  which  the  parameter  is  based 

A,  =  polynomial  coefficients 

i-1  =  order  of  polynomial 
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Fanning  Friction  Factor   (function  of  Re;   1000  *  Re  *  5000) 

•1 
.-4 
■6 
-9 
■12 
-16 


A  =  -1.311518618929068  x  10 

A  =  7.522772995845896  x  10" 

A  =  -1.550552244236629  x  10' 

A,  -  1.730155685156651  x  10" 

A  =  -1.189041200805969  x  10" 


A?  =  -1.590750036125343  x  10 

A0  =  3.138452606608756  x  10 
o 

A9  =  -3.933692644883426  x  10 

A10-  2.838094432805311  x  10 

A  .=  -8.976537263266828  x  10 


-19 


-23 


-27 


■31 


-36 


A,  =  5.329916594794135  x  10 
6 


Colburn-j  Factor  (function  of  Re;   1000  *  Re  <-  5000) 

•1 


A  =  1.644630125351775  x  10 

A  =  -6.796500408006909  x  10 

A3  =  1.238993655152145  x  10 

A,  =  -1.292429783097030  x  10 


-4 


-6 


-9 


A  =  8.543160388163581  x  10 

A,  =  -3.740734237674230  x  10 
b 


-13 


•16 


A  =  1.100442463486636  x  10 

AQ  =  -2.152007282686803  x  10 
o 

A9  =  2.683333339022575  x  10 

A  =  -1.930653368061576  x  10 

A  -  6.099637443396928  x  10 


-19 


-23 


-27 


•31 


-36 


Prandtl   Number    (function  of   temperature) 


A     =  0.7865630968 

A     =-0.2871584338  x  10 

A3  =  0.3666074673  x  10 

A,  =-0.1373281188  x  10 
4 

A  =  0.1569904076  x  10 


-4 


-fc 


-9 


-11 


A^  =  -0.1973361281  x  10 
b 

A  =  0.1041204710  x  10 

A_  =  -0.2047252924  x  10 
o 

A  =  -0.1538501631  x  10 

A  =  0.7960530427  x  10 


-14 


■17 


-21 


-25 


■29 


Viscosity      (function  of   temperature) 


10, 


JUL 


;    (Q)    x    (1.716  x   10        )    x  3600.0 


B.2 


14.822 


where  Q   is   determined  by   the   polynomial  equation 
(as  a   function  of   temperature)    of  coefficients: 
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A  =  -0.4132588012  x  10 

A  -  0,2706898585  x  10 

A  =  -0.1318733093  x  10' 

A,  =  0.7538864121  x  10 
4 

A  =  0.7221455415  x  10 


-1 


-2 


-12 


-12 


A6  = 
A7  = 
A8  = 
\  = 
A10= 


-0.5752187450  x   10 


0.2117993215  x   10 


■0.3862860359  x    10 


0.2782105119  x   10 


0.6763288488  x   10 


■15 


18 


■22 


■26 


■32 


Specific  Heat      (function  of   temperature) 
A     -     0.24913875  x   10° 


A     =   -0.48196782  x   10 

A     =     0.68100588   x   10 
3 


-7 


A4  = 
A5  = 
A6  = 


-0.21358038  x  10 


0.48224330  x  10 


0.47985777  x  10 


-10 


-15 


-18 


Enthalpy   (function  of  temperature) 

2 


A  =  -0.97746770  x 
A  -  0.25257199  x  10 
A3  =  -0.31167200  x  10 


0 


-4 


A,  =  0.30219636  x  10 
4 

A5  =  -0.94272843  x  10" 

A,  =  0.11868992  x  10 ' 
6 

A  =  -0.32615268  x  10 


■11 


•14 


-19 


Relative  Pressure   (function  of  temperature) 


P  ,  -    0.1  Lxp  (D) 

rel 


B.3 


where  D  is  determined  by  the  polynomial  equation 
(as  a  function  of  temperature)  of  coefficients: 


A  -  -0.4043080  x  10 

A  =  0.2759446  x  10" 

A  -  -0.45085958  x  10' 

A.    =  0.53019327   x   10° 
4 

A     -  -0.41468032  x   10" 


•10 


kr    =     0.21127191  x   10 


-13 


A  =  -0.67216409  x  10 

A_  =  0.12101712  x  10 

o 

A     =  -0.94003173   x  10 


-17 


-20 


■25 
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Temperature   (function  of  enthalpy) 

The  temperature  is  computed  in  the  standard  polynomial  form,  as  a 
function  of  x,  where 


A  =  0.23796617  x  10 

A2  =  0.37807000  x  10 

A  -  0.49211414  x  10' 

A,  =  -0.15095696  x  10' 

A  =  0.36226837  x  10 


x  =  log(10.0 
3 


temperature) 


A,  =  -0,45916847  x  10 
6 

A?  =  0.33772184  x  10 

A  =  -0.13466334  x  10 

A  =  0.23704971  x  10 


B.4 


-1 


-2 


-4 


Temperature   (function  of  enthalpy) 
A  =  0.40000290  x  103 
A2  =  0.41827332  x  101 


A  =  -0.64799748  x  10 

A,  =  0.12160439  x  10 
4 

A  -  -0.15484329  x  10 


-4 


-6 


A6  = 
A7  = 
A8  = 
A9  = 


0.86703630  x  10 


-9 


-0.27042766  x  10 


0.48670442  x  10 


■11 


•14 


-0.47242455  x   10 


10 


-  0.19149538  x   10 


-17 


-20 
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/ 

Assumed  Initial  '" 

values  of  m,  £  , 

A  recuperator 

pressure  drops 

I 

/r 


Fixed  cycle 
parameters 


X 


% 


Calculate  thermodynamic  properties  at  cycle 
points  2  4  5,  S  P  ,  I  ,  &  P^ 


Calculate 7 

t3,  i5>  y  a-9 


Recalculate 
recuperator  pressure  drops 


X 


Yes  .. 


^      is  new  P  within 
tolerance  level  of  old 


V 


>-  No 


Recalculate  cycle  parameters 
based  on  new  P 


Calculate  r^   &  HP 


Yes 


I 


Is  HP 
decreas  i  ng? 


•Desired  HP 
Not  Obtainable' 


No 


Yes 


Is  HP  within 
tolerance  level  of 
desired  HP? 


-Yes 


No  ^ 


Is  i~kh  within 
tolerance  level  of 
desired  fkv*  ? 


Calculate  new  m 
from  present 
parameters 


\Prlnt  Outputs/ 


(|top) 
FIGURE  25    PROGRAM  FLOW  CHART 
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'Afrl'Afrmax'Ll'Lmax' 
desl  red  ^Kt    A  &  L 
Increments 


L  -  L 


Yes 


Calculate  outputs  in 
main  subroutine 

Is  HP  within 
tolerance  level  of 
desired  HP? 


Yes 


Is  ru*  within 
tolerance  level  of 
des  I  red  rUvv  ? 


No 


p(lsruv\    decreaslng?)- 


Yes 


No 


yruh  unobtainabley     j(l  s  n^  too  low?) 


Yes 


Increment  L 


I 


No 


Decrease  L. 

■  »  6 


r(ls  L  too  large? 
Yes" 

»L  >  L    \7 

max  / 


No 


Print  outputs 


-* — o. 


Increment  A 


fr 


Yls  A.   >  A-   "TV, 
\L  tl  frmax  / 


Yes  " 


No 


(Sto£) 
FIGURE  26    OUTPUT  III  FLOW  CHART 
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'     bhl'bhmax'bh 
increments  & 
Output  III  Inputs 


bh  -  bhi 


I 


Compute  £>v\  &  r 


Print  ftK*  l~j 


Go  to  Output  I  I  I 


\Prlnt  outputs/ 


\    h    hmax 


Yes 


No 


Increment  b. 


(Stop) 


FIGURE  27    OUTPUT  VIII  FLOW  CHART 
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ZUTPUT  TYPE  THREE 

SPCCIAL  INPUT  PARAMETEFS 

MINIMUM  FRCNTAL  AREA  =  2.3C0  SQ.FT.       LMIN  =  5.350  FT. 

MAXIMUM  FRCNTAL  AREA  =  2.4(3  SC.FT.       LMAX  =20.000  FT. 

FRONTAL  ARE*  INCREMENT  =  0.005  SO. FT.     LINC  =  0.005  FT. 

THERMAL  EFFICIENCY  =(.34f 
OUTPUT  CESIREC  -  CORRESPONDING  LENGTHS 


DESIRED  THERMAL  EFFICIENCY  NOT  OBTAINABLE  WITH  FRCNTAL  AREA  =  2.3C0 

MAXIMUM  THERMAL  EFFICIENCY  =0.337  AT  A  LENGTH  OF   5.350  FT. 

THEPMAl  SPECIFIC  FUEL  TO  AIR      MASS  FLC 

HCRStPUWER      EFFICIENCY  FUEL  CON  RA*no           RATE 

599.92        L.3372  0.397  0.0108          22140. 

TRANSFER  REGENERATOR  AIR  SIDE        GAS  SIDE 

UNITS         CMIN/CMAX  EFFECTIVENESS  REYNOLDS  NC    REYNOLDS  N 

3.00            0.966  0.759  4109.          3751. 

AIR  SIDE        GAS  SIDE  REGENERATOR  RFGENEPATCR   RATIC  GAS/A1 

DELTA  P/P                 DELTA  P/P  VOLUME  WEIGHT      FRONTAL  ARE/1 

0.009           C.C65  12.3050  399.71           1.292 

P?        P3        Pt  P5        T2  T3        T5        T6 

58.89     c3.35     57.18  15.75     873.     1487.  1681.     IC38. 


FRCNTAL  AREA  =  2.3^5  SO. FT.  REQUIRED  LENGTH  =  5.6800  FT. 

THERMAL  SPECIFIC  FUEL  TC  AIR      MASS  FLC 

HCRSEPOWER      EFFICIENCY  FUEL  CON  RATIO           RATE 

5°9.51        C.33c5  0.394  0.0106          22279. 

TRANSFER  REGENERATOR  AIR  SIDE        GAS  SIDE 

UMTS         CHIN/CMAX  EFFECTIVENESS  REYNOLDS  NO    REYNOLDS  N 

3.18            ^.967  C.770  4115.          3772. 

AIR  SIDE        GAS  SIDE  REGENFRATCR  REGENERATCR   RATIC  GAS/AI 

DELTA  P/P       DELTA  &/P  VOLUME  WEIGHT      FRONTAL  AREA 

0.010           C.C6Q  13.0924  425.28           1.292 

P2        P3        P4  P5        T2  T3        T5        T6 

58.89     53.31     57.15  15.81     873.  1497.     1683.     1079. 


FIGURE  28    COMPUTER  OUTPUT 
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APPENDIX  C 


ERROR  ANALYSIS 


It  is  a  basic  fact  that  all  engineering  investigations  are  subject 
to  inaccuracies.   As  in  any  other  study,  a  computer  analysis  is  in- 
herently in  error  due  to:   (1)  the  basic  assumptions  involved  to  facili- 
tate the  problem  solution  (as  neglecting  plate  heat  transfer  resistance 
in  equation  4.7)  and  (2)  the  approximations  made  in  order  to  use  pre- 
viously compiled  engineering  data  (as  the  polynomial  curve  fittings  dis- 
cussed in  Appendix  B) .   A  complete  examination  of  all  the  errors 
introduced  by  these  factors  would  be  not  only  time  consuming,  but  quite 
pointless. 

The  most  significant  error  in  this  computer  analysis  arises  from 
the  use  of  tolerance  limits  on  necessary  output  parameters  and  finite 
increments  on  input  parameters.   This  investigation  necessitates  the 
use  of  two  sets  of  tolerance  levels.   The  first,  Set  I,  involves  a 
tolerance  limit  of  3.3%  (0.3389  ^  r^  1    0.3411)  on  the  desired 
thermal  efficiency  of  0.34  and  a  limit  of  0.2%  (599.0  L    R?      L     601.0) 
on  the  desired  600  horsepower.   Data  for  these  limits  is  calculated  at 
frontal  area  increments  for  the  recuperator  core  of  0.02  sq  ft  and 
length  increments  of  0.01  ft,    Further  investigation,  in  the  area  of 
optimization,  requires  tighter  tolerances  of  1.45%  (0.3395  — Ua^-  0.3405) 
on  the  thermal  efficiency  and  0.09%  (599.5  -  HP  -  600.5)  on  the 
horsepower.   Data  for  these  levels,  designated  as  Set  II,  is  supplied 
at  frontal  area  increments  of  0.005  sq  ft  and  length  increments  of 
0.005  ft.   Data  from  Set  I  is  reproduced  in  Figures  8  through  14  and 
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Tables  III  through  IV;  Set  II  data  is  presented  in  Figures  15  through 
22  and  Tables  V  through  VIII. 

As  there  are  no  performance  characteristics  available  for  the 
majority  of  recuperator  cores  investigated  in  this  study,  no  judgement 
can  be  made  as  to  the  errors  introduced  by  the  computer  simulation.   An 
investigation  of  variations  in  core  parameters  as  presented  in  Table  VIII, 
however,  yields  an  idea  as  to  the  relative  error  introduced  by  the  neces- 
sary use  of  finite  tolerances  and  increments. 

It  should  first  be  noted  in  Table  VIII  that  three  data  sets  are 
presented  for  each  core.   Set  III  involves  the  use  of  tolerance  levels 
as  specified  for  Set  II,  with  reduced  frontal  area  and  length  increment 
magnitudes  of  0.002.   Variations  between  Sets  II  and  III  are  so  small 
that  the  added  computer  calculation  time  introduced  with  the  use  of  Set 
III  over-rides  the  slight  benefits  due  to  the  smaller  increments. 

Variations  in  all  outputs  except  dimensional  and  weight  parameters 
are  of  such  an  order  that  the  error  in  computing  these  quantities  seems 
negligible.   The  most  noticeable  differences  between  output  sets  arise 
in  the  frontal  area  and  length  requirements  and  the  subsequent  weight 
and  volume  quantities.   The  effect  of  tolerance  levels  on  dimensional 
requirements  is  more  noticeable  in  Figures  29  and  30. 

As  may  be  seen,  for  any  given  frontal  area,  the  tighter  tolerance 
of  Set  II  produces  a  core  of  greater  length  than  Set  I.   A  comparison 
of  the  two  plots  results  in  the  fact  that  this  effect  is  more  notice- 
able at  larger  core  spacings.   For  all  cores  tested,  length  variations 
between  sets  decrease  as  frontal  areas  increase  and  less  core  length  is 
required.   It  is  difficult  to  establish  a  trend  as  to  the  frontal  areas 
of  minimum  volume  cores  for  the  two  sets.   Although  Table  VIII  indicates 
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that  Set  II  produces  minimum  volume  cores  of  greater  frontal  area  than 
Set  I,  this  effect  does  not  hold  for  all  cores  investigated.   However, 
all  Set  II  minimum  volume  cores  of  less  frontal  area  than  the  corres- 
ponding Set  I  minimum  volume  core  require  a  much  greater  length  than 
that  for  the  Set  I  core.   Thus,  the  overall  volume  and  weight  trends 
between  cores  for  each  specific  set  seem  to  hold. 

The  most  important  questions  concerning  relative  error  of  the  analy- 
sis have  been  answered.   There  is  evidently  little  error  in  the  perform- 
ance characteristics  output.   Dimensional  comparisons  between  the  two 
sets  result  in  a  maximum  difference  in  length  of  0.1  ft  for  any  given 
frontal  area  and  0.1  sq  ft  in  area  for  any  minimum  volume  core.   Most 
important,  due  to  similar  trends  in  the  two  sets,  the  data  presented 
may  be  assumed  to  provide  specific  patterns  in  plate-fin  counterflow 
recuperator  core  dimensions  and  weight  requirements,  while  the  actual 
magnitudes  presented  may  be  assumed  to  be  within  a  reasonable  range  of 
those  which  might  actually  be  required. 
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